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ABSTRACT 



Massive stars that lose their hydrogen-rich envelope down to a few tenths of a solar mass explode as extended type lib supernovae, 
an intriguing subtype that links the hydrogen-rich type II supernovae with the hydrogen-poor type lb and Ic. The progenitors may be 
very massive single stars that lose their envelope due to their stellar wind, but mass stripping due to interaction with a companion star 
' in a binary system is currently considered to be the dominant formation channel. 

• ' Anticipating the upcoming automated transient surveys, we computed an extensive grid of binary models with the Eggleton binary 

i-C evolution code. We identify the limited range of initial orbital periods and mass ratios required to produce type lib binaiy progenitors. 

The rate we predict from our standard models, which assume conservative mass transfer, is about six times smaller than the current 
rate indicated by observations. It is larger but still comparable to the rate expected from massive single stars. We evaluate exten- 
?^ ' sively the effect of various assumptions such as the adopted accretion efficiency, the binary fraction and distributions for the initial 

, binary parameters. To recover the observed rate we must generously allow for uncertainties and consider low accretion efficiencies in 

^ ■ combination with limited angular momentum loss from the system. 

Motivated by the claims of detection and non-detection of companions for a few lib supernovae, we investigate the properties of the 
secondary star at the moment of explosion. We identify three cases: (1) the companion is predicted to appear as a hot O star in about 
90% of the cases, as a result of mass accretion during its main sequence evolution, (2) the companion becomes an over-luminous 
^ . B star in about 3% of the cases, if mass accretion occurred while crossing the Hertzsprung gap or (3) in systems with very similar 

^SJ ' initial masses the companion will appear as a K supergiant. The second case, which applies to the well-studied case of SN 1993J and 

CO ' possibly to SN 200 lig, is the least common case and requires that the companion very efficiently accretes the transferred material - in 

, contrast to what is required to recover the overall lib rate. We note that relative rates quoted above depend on the assumed efficiency of 

semi-convective mixing: for inefficient semi-convection the presence of blue supergiant companions is expected to be more common, 
occurring in up to about 40% of the cases. 

Our study demonstrates that type lib supernovae have the potential to teach us about the physics of binary interaction and about 
stellar processes such as internal mixing and possibly stellar-wind mass loss. The fast increasing number of type lib detections from 
automated surveys may lead to more solid constraints on these model uncertainties in the near future. 

Key words. Stars: evolution - Binaries: general - Supernovae: general - Supernova: individual: SN 1987K, SN 1993J, SN 1996cb, 
SN 2000H, SN 2001ig, SN 2001gd, SN 2003bg, SN2008ax, Cas A 



1 . Introduction scattered light echo (Krause et al. 2008). The typical light curve 

of a type lib supernova, such as SN 1993 J, is characterized by 
Core collapse supernovae are the bright explosions marking the ^^^y.^ ^he first maximum is associated with shock heat- 
end of the lives of massive stars. Their light curves and spectral jjjg j^e hydrogen-rich envelope, resembUng a type II super- 
signatures come in a large variety of types and yield information ^^^^ ^he second maximum is caused by the radioactive decay 
about the structure and chemical composition of the progenitor nj^kel (Benson et al. 1994). These characteristics can be ex- 
star and It surroundings. Type II supernovae - characterised by pj^jngj assuming that the progenitor star had an extended low- 
strong hydrogen lines - are associated with massive stars that ^ass hydrogen envelope at the time of explosion, with between 
are still surrounded by their hydrogen-rich envelope at the time q.I and 0.5 of hydrogen (Podsiadlowski et al. 1993; Woosley 
of explosion, whereas type lb and Ic supernovae - in which no ^ 1994. Elmhamdi et al. 2006). Progenitors with smaller hy- 
signature of hydrogen is found - are thought to result from mas- jrogen envelope masses are compact, but can be classified as 
sive stars that have lost their entire hydrogen-rich envelope. ^ypg jj^ supernovae if hydrogen lines are detected shordy after 
Type lib supernovae constitute an intriguing intermediate the explosion, e.g. Chevalier & Soderberg (2010). In this paper 
case. Initially they show clear evidence for hydrogen, but later we focus on the progenitors of extended type lib supernovae. 
the hydrogen lines become weak or absent in the spectra. Two Massive single stars with masses of about 30 Mg or higher 
famous examples of this subtype are SN 1993J (Ripero et al. can lose their envelope as a result of a stellar wind. Alternatively, 
1 993) and Cassiopeia A, which was recendy classified using the massive stars can be stripped of their envelope due to interaction 

with a companion star in a binary system. Quickly after the dis- 

* Hubble fellow. CO very of SN 1993 J many authors realized that the single star 
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scenaiio requires precise fine tuning of the initial mass of the 
star in order to have a very low-mass envelope left at the time of 
the explosion (Podsiadlowski et al. 1993; Woosley et al. 1994, 
and references therein). According to these authors, mass strip- 
ping from an evolved red supergiant in a binary system naturally 
leads to a hydrogen envelope containing a few times 0. 1 Mq at 
the moment of explosion. More compact binaries, where the pri- 
mary star is stripped in a less evolved stage, result in smaller 
hydrogen envelopes and may lead to type lib supernovae of the 
compact category (e.g. Yoon et al. 2010) or type lb supemovae. 

Nomoto et al. (1993) proposes a third scenario which in- 
volves a common envelope phase and finally a merger of the 
core of the primary star and its companion. Energy to eject the 
envelope is extracted from the orbital energy, but a small layer of 
hydrogen may remain on the surface of the merger product until 
the moment of explosion as a supernova type lib. This scenario 
was recently considered by Young et al. (2006) for Cassiopeia 
A, in which no evidence for a companion star was found, de- 
spite numerous attempts. 

A decade after the detection of 1993 J, the interest for type 
lib supernova revived when Maund et al. (2004) found evidence 
for the presence of a blue supergiant companion at the loca- 
tion where 1993J had just faded away. These findings nicely 
fitted the predictions for the companion star in the models by 
Podsiadlowski et al. (1992, 1993). StancHffe & Eldridge (2009) 
performed the first systematic study of binary progenitor mod- 
els that matched the properties of the progenitor of 1993 J and 
its companion. They computed a grid of binary models vary- 
ing the initial primary masses and initial orbital periods. Even 
though they find a suitable progenitor model, they emphasize 
that it proves to be extremely difficult to explain the properties 
of the blue supergiant companion star. 

At the time of writing about 69 supernovae have been classi- 
fied as lib supernova'. This number will increase quickly thanks 
to upcoming systematic wide-field surveys of transient events, 
such as the Palomar Transient Factory (e.g. Arcavi et al. 2010), 
Pan-STARRS (Kaiser et al. 2002) and the Large Synoptic Survey 
Telescope (Ivezic et al. 2008). About 9 have been studied well 
enough such that constraints can be derived on the progenitors 
and in some cases on possible companion stars. 

These observational developments and the questions raised 
by Stancliff'e & Eldridge (2009) motivated us to further investi- 
gate binary progenitor models for type lib supernova. In con- 
trast with previous work we will not just discuss the case of 
SN 1993J. Instead we discuss the properties of the progenitor 
stars and their companions in general. We extend the work of 
Stancliffe & Eldridge (2009) by exploring the parameter space of 
binary progenitor models in detail, in particular the dependence 
on the initial mass ratio. Furthermore we discuss the properties 
of the companion star at the time of explosion. 

In section 2 we review observed type lib supemovae and the 
observed rates. In section 3 we outline the assumptions in our 
binary evolution code. In section 4 we briefly discuss single-star 
progenitor models. Section 5 we discuss our progenitor models 
and in section 6 the rate of lib supernovae predicted from our 
models. Section 7 contains a discussion, conclusion and outlook. 

2. Observed type lib Supernovae 

In recent years about 69 supemovae of type lib have been 
detected'. These observations enabled the first estimates for the 
rate of type lib supernovae. The determination of the rate is not 

' http://www.cfa.liarvard.edu/iau/lists/Supemovae.html 



straightforward due to observational biases and selection effects. 
For example, if the supernova is observed too late it will be clas- 
sified as a type lb supernova instead of type lib (see also Maurer 
et al. 2010). Van den Bergh et al. (2005) and Li et al. (2007) 
derive that 3.2 + 1.0% and 1.5 + 1.5% of all core coflapse SNe 
(CCSNe) are of type lib. They based their estimates on the dis- 
coveries by the Lick Observatory Supemova Search (LOSS), re- 
spectively using a 140 and 30 Mpc distance-limited sample. Van 
den Bergh et al. (2005) also specified that the distribution of SNe 
lib and SNe II does not depend on the morphological type of the 
host galaxy. Smartt et al. (2009); Smartt (2009) evaluated dif- 
ferent determinations of the rate and the effect of observational 
biases. They estimate the rate to be 5.4 + 2.7% based on ob- 
servations covering 10.5 years within 28 Mpc. However, there 
are only 5 lib supernovae in their sample, of which at least two 
are compact lib supemovae. Preliminary results of the Palomar 
Transient Factory survey point towards a rate of 3.6 + 2.5% in 
giant hosts and a larger fraction in dwarf galaxies, 20 + 1 1 % 
(Arcavi etal. 2010). 

While this study was near completion Smith et al. (2010) 
presented the results of a homogeneous volume limited survey 
within 60 Mpc which includes 80 CCSNe. They find a signif- 
icantly higher fraction of type lib supernova, 10.6;^3 j%, with 
respect to previous studies. They attribute the difference to their 
more complete photometric and spectroscopic follow-up obser- 
vations used to classify type II supernovae into different sub- 
types. It is unclear how many of the type lib supemovae quoted 
in these studies are extended instead of compact. 

For the purpose of comparison with our model predictions 
we will use the rate quoted by Smartt et al. (2009) after taking 
out the two compact Type lib, resulting in a fraction of extended 
type lib of about 3% with respect to all core collapse supernovae. 
We emphasize that this number is still very uncertain. More reli- 
able determinations of the rate of extended type lib supemovae 
are expected in the near future from the current and upcoming 
automated transient surveys." 

In the follow paragraphs we discuss some observed super- 
novae type lib and discuss their observed characteristics. 

2.1. Individual type lib supernovae 

2.1.1. SN 1987K 

SN 1987K was the first supernova observed with the characteris- 
tics of a type lib supernova, namely the transition from a super- 
nova type II to type lb (Filippenko 1988), but it was only later 
that it was defined a SN of type lib. Filippenko (1988) already 
suggested that these characteristics could be due to mass loss of 
a massive star, with an initial mass between 20 and 25 M© or a 
combination of mass transfer and winds of a less massive star, 
between 8 and 20 M©. 

2.1.2. SN 1993J 

SN 1993J was the first supernova classified as type lib. This 
supernova has been well studied because the progenitor was 
detected and recognized as a star with a spectral class KOIa 
(Filippenko et al. 1993). Its binary companion has been observed 
recently as an early B-supergiant (with best estimate a B2Ia star, 
Maund & Smartt 2009). The light curve of the supernova can 
be explained if the star had an amount of hydrogen in the enve- 
lope of about 0.1-0.5 M© at time of explosion (Woosley et al. 
1994; Filippenko et al. 1993). Maund et al. (2004) determined 
the luminosity and effective temperature of the progenitor of the 
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supernova and its companion, namely log L/Lq = 5.1 + 0.3 and 
log Tetf/K = 3. 63 +0.05 for the progenitor and log L/Lq = 5+0.3 
and log Teff/K = 4.3 +0.1 for the companion. The ejecta mass 
was determined, using the width of the second peak of the light 
curve, to be about 4 Mq (Shigeyama et al. 1994). Radio and X- 
ray observations showed evidence for a mass loss rate of about 
4x10"^ Mq yr"' at the time of explosion (Fransson et al. 1996). 

2.1.3. SN 1996cb 

Supernova 1996cb was also determined to be of type lib (Qiu 
et al. 1999). The light curve of this supernova was not observed 
until several days after the explosion. Instead of the two peaks 
that are typical for type lib light curves, the light curve con- 
sists of a short-term plateau phase similar to a type II supernova. 
These differences with SN 1993 J arise from a more massive hy- 
drogen envelope of SN 1996cb (Qiu et al. 1999). 

2.1.4. SN2000H 

Benetti et al. (2000) classified SN 2000H a type lib super- 
nova because of the hydrogen lines observed in the spectrum. 
However, these hydrogen lines were not obvious and the light 
curve showed more resemblance with a SN lb (Branch et al. 
2002). Therefore, this supernova has been classified by other 
authors as a type lb supernova (Branch et al. 2002; Elmhamdi 
et al. 2006). The amount of hydrogen was estimated to be 0.08 
Mo (Elmhamdi et al. 2006). 

2.1.5. SN2001gd 

The observations of radio light curves give an estimate of the 
mass loss rate of the progenitor system of the supernova. The 
radio light curve of the type lib SN 2001gd indicates a mass loss 
rate of about 2-12 x 10"^ Moyr"' (Stockdale et al. 2003; Perez- 
Torres et al. 2005), which is in between the rates for typical type 
II and type lb SNe. 

2.1.6. SN2001ig 

Another supernova of type lib, SN 200 1 ig, was observed in 200 1 
in NGC 7424 and shows a spectral evolution similar to that of 
SN 1993J (Ryder et al. 2006). Evidence was found for a star of 
spectral type late-B through late-F at the location of SN 200 lig, 
a possible companion of the progenitor of SN 200 lig (Ryder 
et al. 2006). 

2.1.7. SN2003bg 

SN 2003bg evolved from a type Ic supernova to a hydrogen-rich 
type lib, to a hydrogen-poor type Ibc (Soderberg et al. 2006). 
It was observed as a broad-lined type lib supernova and pro- 
claimed to be 'the first type lib hypemova' (Mazzali et al. 2009; 
Hamuy et al. 2009). The broadness of the lines indicates a high 
progenitor mass (Hamuy et al. 2009). The light curve and spec- 
tral evolution indicate the presence of a thin layer of hydrogen 
at time of explosion, x 0.05 Mq (Mazzali et al. 2009). The ve- 
locity of the ejecta resembles more closely the velocity of SNe 
type lb than type II (Soderberg et al. 2006). This implies a com- 
pact Wolf-rayet progenitor, with a progenitor mass between 20 
and 25 M©. Soderberg et al. (2006) conclude that this event is 
an intermediate case between SNe type lib and type lb. 



2.1.8. SN 2008ax 

The light curve of SN 2008ax shows some diff'erences with the 
light curve of SN 1993 J, namely the lack of the first peak and it 
has slightly bluer colors (Pastorello et al. 2008). These features 
can be explained by a less massive hydrogen envelope at time 
of explosion, less than a few x 0. 1 M©, in comparison with the 
progenitor of SN 1993J (Crockett et al. 2008). 

2.1.9. Cas A 

Cas A is the supernova remnant of a star that exploded about 
350 years ago (Thorstensen et al. 2001). A light echo from this 
explosion (Krause et al. 2008) shows evidence that it was a su- 
pernova of type lib. Direct methods to determine the mass of the 
progenitor star are difficult, but the ejecta mass was calculated to 
be 2-4 Mq and the remnant would be expected to be a neutron 
star with a mass between 1.5 and 2.2 M© (Young et al. 2006). 
This sets the mass of the star at time of its explosion at about 4-6 
Mq. There is no direct evidence as to wether this supernova was 
of the compact or extended type lib, but the possibilty that the 
progenitor was a red supergiant is left open. For this supernova 
single and binary progenitor models were calculated (Young 
et al. 2006). The single star models indicated fine-tuning of the 
stellar wind is necessary to evolve to the specific characteristics 
of the supernova remnant Cas A (Young et al. 2006). Besides, 
there is evidence that the progenitor could only have had a 
very short-lived Wolf-rayet phase, which is difficult to explain 
with single stars (Schure et al. 2008; Van Veelen et al. 2009). 
There is no evidence for a companion star. Therefore a common 
envelope scenario was proposed, in which the two stars merge 
into a single star before explosion. Observations show tentative 
evidence for this scenario (Krause et al. 2008), such as the 
asymmetric distribution of the quasi-stationary flocculi near Cas 
A, which could arise from the loss of a common envelope. 

The observations put constraints on the general properties 
of a SN type lib: the explosion of a supergiant with a hydrogen 
envelope mass between 0.1 and 0.5 Mq. We consider the lower 
limit of the mass of the hydrogen envelope to be 0. 1 M© rather 
than 0.01 Mq, the lower limit proposed by Chevalier & 
Soderberg (2010). The explosion of a star with a hydrogen 
envelope mass smaller than 0. 1 Mq will exhibit hydrogen lines 
in its spectrum, but only in the early phases of the supernova. 
In addition the light curve resembles the typical light curve of 
a SN type lb. Therefore such a supernova will more likely be 
defined as a supernova type lb or a transitional type between lib 
and lb (see examples above, e.g. SN 2003bg and SN 2000H). 
Elmhamdi et al. (2006) places the upper limit of the mass of the 
hydrogen envelope of the progenitor of a type lb supernova at 
0.1 Mo. 

A binary progenitor is confirmed or considered likely in 
some cases, and the secondary has been detected as a blue su- 
pergiant in possibly two cases. However, no general constraints 
on the secondary can be set. 

3. Stellar evolution calculations 

We use a version of the binary evolution code STARS origi- 
nally developed by Eggleton (1971) and later updated and de- 
scribed by various authors (e.g. Pols et al. 1995; Eggleton 2006; 
Glebbeek et al. 2008). The code is fully implicit and solves the 
equations for the structure and composition of the star simultane- 
ously. It employs an adaptive non-Lagrangian mesh that places 
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mesh points in regions of the star where higher resolution is re- 
quired. This allows us to evolve stars with a reasonable accuracy 
using as few as 200 mesh points. The code therefore is fast and 
suitable to compute the large numbers of models needed to in- 
vestigate wide initial parameter space of binary systems (e.g. De 
Mink et al. 2007). 

We use nuclear reaction rates from Caughlan et al. (1985) 
and Caughlan & Fowler (1988) and opacities from Rogers & 
Iglesias (1992) and Alexander & Ferguson (1994). The as- 
sumed heavy-element composition is scaled to solar abundances 
(Anders & Grevesse 1989). 

Convection is implemented using a diffusion approxima- 
tion (Eggleton 1972) of the mixing-length theory (Bohm-Vi tense 
1958), assuming a mixing length of 2.0 pressure scale heights. 
We use the Schwarzschild criterion to determine the boundaries 
of the convective regions. Convective overshooting is taken into 
account using the prescription of Schroder et al. (1997) with 
an overshooting parameter of 6ov - 0.12, which was calibrated 
against accurate stellar data eclipsing binaries (Pols et al. 1997). 
In terms of the pressure scale height, as the overshooting param- 
eter is commonly defined in other stellar evolution codes, this 
value approximately compares to ffov ~ 0.25. 

Mass loss in the form of a stellar wind is taken into account 
adopting the prescription by de Jager et al. (1988). Although for 
O and B stars this prescription has been superseded by more re- 
cent mass-loss determinations, this is not the case for the red 
supergiant region of the H-R diagram where the only signifi- 
cant mass loss in our binary models occurs. For computational 
reasons we ignore stellar-wind mass loss from the less massive 
companion star and we ignore any possible accretion from the 
stellar wind of the primary. 

To compute the evolution of interacting binaries we evolve 
the two stars quasi-simultaneously. First we follow the primary 
for several steps and store the changes of the masses of both 
stars and of the orbit. Afterwards the secondary is evolved ap- 
plying these mass changes until its age reaches that of the pri- 
mary. When the primary star expands beyond its Roche lobe we 
compute the mass-transfer rate as a function of the difference in 
potential between the Roche-lobe surface and the stellar surface. 
The mass flux of each mesh point beyond the Roche lobe is given 
by 



(1) 

as in De Mink et al. (2007), where m and r denote the mass and 
radius coordinate of the mesh point and 0s the difference in po- 
tential with respect to the Roche-lobe surface. C is a proportion- 
ality constant which we set to 10"^ for numerical convenience. 
The mass transfer rate is then given by the integral of Eq. (1) 
over all mesh points outside the Roche radius. As a result of the 
implementation of mass transfer in our code we implicitly as- 
sume that the entropy and composition of the accreted material 
is equal to the entropy and composition of material at the surface 
of the accreting star (e.g. Pols 1994). Tidal effects are not taken 
into account in this work. 

To investigate the effects of non-conservative mass trans- 
fer, i.e. mass and angular momentum loss from the system dur- 
ing Roche-lobe overflow, we assume that a constant fraction 
/3 of the transferred mass is accreted by the companion. We 
adopt different values: P - I (conservative mass transfer) and 
P - 0.75,0.5,0.25 (non-conservative mass transfer). The spe- 
cific angular momentum of the mass lost from the system is as- 
sumed to be the specific angular momentum of the orbit of the 
accreting star. 



Table 1. Properties of single stellar models for different initial masses 
Mi. 



Mi 

(Me) 


Mf 
(Mo) 


Mh 
(Mo) 




Tcff 
(lO'K) 


SN type 


32 


13.94 


1.014 


5.56 


3.6 


IIL 


32.5 


12.89 


0.505 


5.57 


3.8 


lib 


33 


11.81 


0.080 


5.58 


4.3 


lb 



Notes. The final mass Mf and the amount of hydrogen in the envelope 
Mh are given at the time of explosion, together with the luminosity 
log(L/Lo) and effective temperature Tefi at the onset of central carbon 
burning. 



Using this code, we calculated a grid of binary systems for an 
initial primary mass of 15 Mq. We varied the initial orbital period 
between 800 and 2100 days in steps of 100 days and we varied 
the initial secondary mass between about 10 and 15 M© in steps 
of 0.1 Mq. For systems with similar masses we increased the 
resolution by varying the secondary mass in steps of 0.01 M©. 
In addition we computed several systems with different primary 
masses. 

The stars are evolved until the onset of carbon burning. The 
star will explode shortly after this point, about 2.6 x 10^ years 
for a star of 15 Mq (El Eid et al. 2004). To determine the masses 
of the stars and in particular of their envelopes at the time of 
explosion we extrapolate the mass-loss rate by stellar winds and 
Roche-lobe overflow to determine the masses at time of explo- 
sion. When both stars fill their Roche lobe at the same time and 
a contact binary is formed, we end our simulation. In this case 
the binary system will probably evolve into a common envelope. 

We note that whether or not a contact system forms is sen- 
sitive to some of our assumptions regarding mass transfer The 
chosen value of C in eq. (1) leads to a maximum mass-transfer 
rate that is lower than the self-regulated rate on the thermal 
timescale of the donor star If we choose a larger value of C 
the mass-transfer rate increases, which affects the response of 
the secondary star Since mass transfer is faster than the ther- 
mal timescale of the secondary, more rapid accretion results in 
a stronger radius expansion of the secondary and a higher likeli- 
hood of forming a contact system. On the other hand, the trans- 
ferred material comes from the surface of a red supergiant which 
has a much smaller specific entropy than the hot surface of the 
secondary. Although the gas may undergo additional heating 
during accretion, this is probably not sufficient to make its en- 
tropy equal to the surface entropy of the secondary as we im- 
plicitly assume. Taking this into account properly, which is very 
difficult, would likely result in less radius expansion of the sec- 
ondary. These two simplifying assumptions thus have opposite 
effects, making the true boundary between systems that do and 
do not come into contact hard to predict from our models (see 
Section 5.2 and Section 7). 



4. Single star progenitors 

At the time of explosion the progenitor of a type lib supernova 
is surrounded by a low-mass hydrogen envelope. In massive sin- 
gle stars the envelope can be removed by the stellar wind. Heger 
et al. (2003) find that at solar metallicity the range for type lib 
and IIL supernovae combined ranges from roughly 26-34 Mq. 
Type IIL supemovae show a linear decline in the light curve, 
but without transition to a type lb, which requires a hydrogen- 
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Mass transfer rate 






Wind mass loss rate 

1 1 1 1 





12 10 8 6 4 2 

Envelope mass {M,J 

Fig. 1. Wind mass loss and mass transfer rate of as a function of the 
remaining envelope mass for the most massive star of binary system 
with initially 15+14.35 Mq, an initial orbital period of 1500 days. For 
details see Section 5. 



rich envelope of less than about 2 M© at the time of explosion. 
Eldridge & Tout (2004) map single-star progenitors of differ- 
ent types of supernovae as a function of metallicity. At solar 
metallicity they find that IIL and lib supernova progenitors com- 
bined should have initial masses between approximately 25 and 
30 Mq. Georgy et al. (2009) discuss single star SN progenitors 
from rotating models at different metallic! ties. Although they do 
not discuss type lib or IIL in particular, they find that the tran- 
sition from type II to lb should occur around 25 Mq. Perez- 
Rendon et al. (2009) investigate possible single star progenitors 
for Cas A and need a progenitor of approximately 30 M©. In par- 
ticular, a single star scenario has been proposed for hypemova 
SN 2003bg to explain the high mass. On the other hand a single 
star scenario cannot explain the characteristics of SN 1993J (see 
Sec. 2.1.2). Consequently, the single star scenario cannot explain 
all observed type lib supernovae. 

The determination of the rate from single stars is sensitive to 
uncertainties in the stellar wind mass loss rates. Especially with 
respect to stellar winds from massive red giants our understand- 
ing is limited. Yoon & Cantiello (2010) speculate that a super- 
wind driven by pulsational instabilities may drive a strong mass 
loss, bringing the minimum mass for type lib supernovae down 
to about 20 Mq. 

Single stars with a mass above ~25 are believed to pro- 
duce only faint supernovae (Fryer 1999). Consequently, these 
type lib SNe will appear different than type lib SNe formed 
by binary stars. Nevertheless, in the correct mass range, single 
stars can explode as type lib SNe and therefore it is reasonable 
to compare the expected rate from binary and single-star pro- 
genitors. To be able to do this comparison we computed single 
stellar models with the same input physics as the binary mod- 
els we discuss in the next section. We find that the initial mass 
range for single stars resulting in type lib progenitors should be 
within 32.5-33 M©. In Table 1 we list the final mass and and the 
amount of hydrogen left in the envelope near this range of initial 
masses. Assuming a Kroupa initial mass function (Kroupa 2001) 
we estimate that about 0.3% of all single stars more massive than 
8 solar masses are within this mass range and will end their lives 
as a type lib supernova. 



5. Binary star progenitors 

In this section we discuss binary progenitor models for type lib 
supernova. We compute the evolution of binary models with dif- 
ferent initial orbital periods and different initial mass ratios. We 
adopt an initial primary mass of 15 Mq in agreement with the 
progenitor model proposed for 1993J (Maund et al. 2004). We 
assume that type lib supemovae result from massive stars that 
undergo core collapse with an envelope which contains between 
0.1 and 0.5 Mq of hydrogen. This criterion is based on obser- 
vations, previous models (Podsiadlowski et al. 1993; Woosley 
et al. 1994; Elmhamdi et al. 2006) and some test models which 
proved that a hydrogen mass less than 0. 1 Mq gives rise to a 
compact rather than an extended progenitor 

As an example we discuss a system with an initial orbital 
period of 1500 days and initial masses of 15 and 14.35 Mq for 
the primary and secondary star respectively. The initially most 
massive star evolves faster and experiences significant mass loss 
in the form of a stellar wind when it ascends the giant branch and 
during central helium burning. After about 13.1 Myr, when the 
helium mass fraction has dropped below 0.5 in the center, it fills 
it Roche lobe (late case B mass transfer, Kippenhahn & Weigert 
1967). At this moment it has already lost more than 1 Mq and 
has become less massive than its companion. The reversal of 
the mass ratio before the onset of Roche-lobe overflow helps to 
stabilize the mass transfer 

In Figure 1 we depict the mass-transfer rate as a function of 
the remaining envelope mass. We find that the mass transfer ini- 
tially takes place on a timescale equal to the thermal timescale 
of the primary star The maximum mass-transfer rate during this 
phase is 6 X 10""* Moyr"'. Although this phase lasts only about 
0.05 Myr, about 4.5 Mq is transferred. After this phase the star 
keeps filling its Roche lobe and mass transfer continues on the 
nuclear timescale, at a rate comparable to the mass-loss rate in 
the form of a stellar wind, about 3 x 10"^ Moyr"'. During this 
phase the primary expands on its nuclear timescale while it is 
burning helium in its center This phase lasts about 0.8 Myr and 
about 2.5 Mq is transfeiTed. After central helium exhaustion the 
star expands again on its thermal timescale and a second max- 
imum in the mass-transfer rate occurs. Finally, at the onset of 
carbon burning the star expands again resulting in a third peak 
in the mass-transfer rate, see Fig. 1. We follow the evolution of 
the system up to this point, when the mass of hydrogen in the 
envelope has decreased to 0.46 Mq. Extrapolating the mass-loss 
rate we find that the amount of hydrogen in the envelope at the 
time of explosion will be about 0.37 Mq. Therefore we expect 
that the primary star explodes as a type lib supernova. 

In wider systems mass transfer starts in a later phase of the 
evolution of the primary star Because the primary stars in these 
systems are more evolved, there is less time available to reduce 
the mass of the envelope before the explosion. In addition, stel- 
lar winds had more time to reduce the mass of the primary star 
before the onset of Roche-lobe overflow. Reversal of the mass 
ratio stabilizes the process of mass transfer This results in a 
lower mass-transfer rate during the first phase of mass transfer 
Vice versa we find that stars in binary systems with lower ini- 
tial orbital periods remain with smaller envelope masses at the 
time of explosion. We give details of all our computed models in 
Tables A.1-A.4 in the Appendix. 

5.1. Properties of the companion 

In general the companion star is relatively unevolved, i.e. still on 
the main sequence, at the onset of mass transfer However, when 
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Table 2. Properties of binary stellar models for an initial primary mass of 15 Mq, conservative mass transfer and with variation of the initial 
secondary mass Mai and initial orbital period FoA- 





M2i 

(Mo) 


(Mo) 


M2f 

(Mo) 


f orb 

(days) 


Mhi 
(Mo) 


log 




(lO^^K) 


log 


[Lei 


(lO'K) 


#1 


14 


5.93 


20.92 


1500 


0.354 


5.07 


3.48 


4.95 


31.4 


#2 


14.35 


5.98 


21.22 


1500 


0.370 


5.07 


3.47 


5.06 


17.9 


#3 


14.55 


6.01 


21.40 


1500 


0.390 


5.07 


3.45 


4.85 


3.82 


#4 


14.95 


6.02 


21.79 


1400 


0.401 


5.07 


3.51 


5.07 


3.63 



Notes. The final mass of the primary and secondary, Mfi and Mf2, and the amount of hydrogen in the envelope Mh are given at the time of 
explosion, together with the luminosity log(L/Lo) and elfective temperature Teff of the primary and secondary at the onset of central carbon 
burning in the primary (except for model #4, where L2 and Teff j refer to the last computed model of the secondary). 



Table 3. Expected broad-band magnitudes of the three dilferent possible evolutionary paths of the companion star, with an intial mass Ma;. 





Mji 
(Mo) 


M(U) 


M(B) 


M(V) 


M(R) 


M(I) 


M(J) 


M(H) 


M(K) 


#1 


14 


-5.89 


-4.99 


-4.65 


-4.66 


-4.37 


-4.02 


-3.91 


-3.80 


#2 


14.35 


-7.12 


-6.49 


-6.22 


-6.27 


-6.03 


-5.77 


-5.68 


-5.65 


#3 


14.55 


-2.45 


-4.66 


-6.14 


-7.45 


-8.15 


-9.07 


-9.78 


-10.2 



the initial mass ratio is close to one, the evolutionary timescales 
of the primary and secondary star are comparable and the com- 
panion can be more evolved. The response of the companion to 
mass accretion depends on its evolutionary stage. We distinguish 
three different cases: (1) accretion starts while the companion is 
on the main sequence, (2) accretion starts while the companion 
is crossing the Hertzsprung gap and (3) accretion starts while 
the companion is a giant. The evolution of the stars in these 
three cases is illustrated in Figure 2, where we give the evolu- 
tionary tracks of both stars in the Hertzsprung-Russell diagram 
and a Kippenhahn diagram illustrating the changes in the inter- 
nal structure as a result of mass transfer Table 2 list several prop- 
erties of the binary models described here. 

The first case, accretion during the main sequence, is de- 
picted in the left panels of Figure 2. The accreting star responds 
to the increase in mass by adapting its internal structure. The 
size of its convective core increases and fresh hydrogen is mixed 
towards the center, effectively rejuvenating the star After accre- 
tion the properties of the star are similar to the properties of a 
younger single star of the same mass. The star becomes brighter 
but remains hot, appearing as an O star, see also Table 2 (model 
#1). 

The central panels of Figure 2 show an example of the sec- 
ond case, accretion during the Hertzsprung gap (model #2 in 
Table 2). In this phase nuclear burning takes place in a shell 
around the core. This prevents the star from adapting its inter- 
nal structure to that of a single star of its new increased total 
mass. Having a core mass which is too small compared to the 
core mass of a normal single star, the star appears as an over- 
luminous B supergiant. This type of progenitor model has been 
proposed to explain the properties of the blue companion of SN 
1993J. 

For systems with very similar initial masses, accretion takes 
place while the secondary resides on the giant branch, see the 
right panels of Fig 2. In this third case we find that the sec- 
ondary will appear as a K supergiant at the moment the primary 
explodes (model #3 in Table 2). If the initial mass ratio is even 
closer to one we find that the evolution of the secondary is accel- 
erated enough to catch up with and overtake the evolution of the 
primary. In this case the secondary explodes first, as a normal 



type II supernova, while the primary explodes afterwards as a 
type lib supernova. The time difference between the explosions 
is 8000 years for model #4 in Table 2, but can be up to 10^ years 
for very close mass ratios (see Table A. 3 in the Appendix). 

While the three companions of the SN lib progenitors in 
these examples are of similar luminosity, it is important to 
consider how easily they might be observed in pre- and post- 
explosion images. For example the O star has a high surface 
temperature and so most of its emission is in the ultraviolet. In 
Table 3 we list the expected broad-band magnitudes of our ex- 
ample progenitors using the methods outlined in Eldridge et al. 
(2007) and Eldridge & Stanway (2009) to calculate the colours. 
We see that while all three progenitors are of a similar bolomet- 
ric luminosity, the B and K supergiants output more of their light 
in the optical bands, such as V and I which are typical of those 
most commonly available in pre-explosion imaging (Smartt et al. 
2009), making them easier to identify in pre- and post-explosion 
images. 

5.2. Parameter space 

Figure 3 depicts the range of initial mass ratios and initial orbital 
periods of binary systems in which the primary star is expected 
to explode as a type lib supernova. In systems with initial or- 
bital periods larger than about 1600-1800 days, depending on 
the mass ratio, the primary star has more than 0.5 Mq of hydro- 
gen left in its envelope at the time of explosion. We assume that 
the supernova would be classified as type II. We therefore do not 
find any SNe type lib progenitors which have undergone Case 
C mass transfer (Lauterborn 1970), since these all end up with 
hydrogen masses greater than 0.5 M©. In systems with initial 
orbital periods smaller than about 1000-1300 days the envelope 
mass left at the time of explosion is less than 0.1 and we 
assume that the primary explodes as a type lb supernova. In sys- 
tems with M2/M1 < 0.7-0.8, depending on the orbital period, 
the mass-transfer rate is so high that the stars come into contact. 
These systems are expected to experience a common envelope 
phase. As we have discussed in section 3, whether or not a bi- 
nary evolves into contact during this phase is sensitive to some 
of our model assumptions. Therefore the critical mass ratio sep- 
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Fig. 2. Response of the accreting star in three different cases (1) left panels: accretion during main sequence, (2) central panels: accretion during 
Hertzsprung gap and (3) right panels: accretion during giant branch. In the top row we show the evolutionary tracks of both stars in the Hertzspmng- 
Russell diagram. In the bottom row we illustrate the evolution of the inteiiral structure of the accreting star as a function of time around the moment 
of accretion. The total mass (black line), the core mass (black line) and the mass coordinates of the region in which nuclear burning takes place 
(red line) are plotted versus time. Grey areas indicate convective regions. In all models we assumed a primary mass of 15 Mg, an initial orbital 
orbital period of 1500 days and a secondary masses of 14, 14.35 and 14.55 Mq to illustrate the three different cases. These diagrams correspond 
to models #1, #2 and #3 in Table 2. 



arating contact from non-contact systems is uncertain and we 
regard the location of this boundary in Fig. 3 as indicative only. 

The borders between type lib / II and between type lb / lib 
run diagonally across this diagram, i.e. the critical orbital pe- 
riod dividing different supernova types increases with mass ratio. 
This is caused by the fact that systems with more extreme initial 
mass ratios exhibit a higher mass-transfer rate in the initial phase 
of mass transfer For a given initial orbital period this results in 
lower envelope masses at the time of explosion in systems with 
more extreme mass ratios. 

The border between type lb and type lib shows a horizontal 
step near mass ratios of about 0.89. This feature is related to the 
evolutionary stage of the primary star at the onset of mass trans- 
fer. For systems with orbital periods smaller than about 1100- 
1200 days, the primary star fills its Roche lobe relatively early, 
while it is still in the Hertzsprung gap just before ascending the 
giant branch. The effects of stellar winds, which enlarge or even 
reverse the mass ratio and stabilize the mass transfer, are still 
limited. In addition, the rapid expansion of the primary on its 
thermal timescale results in a high mass-transfer rate. We find 
that within this period range only in systems with very similar 
initial masses, MnjMx > 0.89, enough hydrogen can be retained 
on the surface of the primary star until the time of explosion to 
result in a type lib supernova. 

In Figure 3 we also indicate the three different cases dis- 
tinguishing the properties of the secondary star at the mo- 
ment of explosion, discussed in Sect. 5.1. The information from 
our model grid on which this figure is based can be found in 



Tables A.1-A.3 in the Appendix. In the most common case 1, 
the companion is still on the main sequence at the onset of ac- 
cretion and it will appear as an O star at the time of explosion. 
Case 2 is the region where the companion accretes while it is 
on the Hertzsprung gap and will appear as a B supergiant. The 
range of mass ratios for this case is very limited. This is a direct 
result of the short time spent by the star in the Hertzsprung gap, 
about 0.5% of time that it spends on the main sequence. Case 3 
indicates the region where the companion will evolve to a K su- 
pergiant. This case is more likely to occur for wider systems in 
which the primary star fills its Roche lobe in a later stage of 
helium burning. In Figure 4 we show an enlarged region of the 
parameter space for initial mass ratios near 1. Here we indicate 
case 3b in which the secondary explodes before the primary. This 
only occurs if the stars are initially very close in mass, to within 
0.5%. 



5.3. Non-conservative mass transfer 

So far we have assumed that the mass transferred during Roche- 
lobe overflow is efficiently accreted by the secondary. This as- 
sumption may not be valid, given the high mass-transfer rates 
reached and the fact that the secondary star may quickly spin up 
to break-up rotation after accreting only a fraction of the trans- 
ferred mass (e.g. Packet 1981). In this section we investigate the 
effect of variations in the accretion efficiency /3, i.e. the fraction 
of material accreted by the secondary with respect to the amount 
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Fig. 3. Ranges of initial mass ratios and orbital periods of binaries for which the primary star is expected to explode as a type lib of the extende type 
supernova. We assumed a primary mass Mi = 15 Mq, a metallicity Z = 0.02 and conservative mass transfer. We assume that type lb supernovae 
result from stars with less than 0.1 Mq of hydrogen in their envelope at the time of explosion, type lib supernovae from stars with 0.1-0.5 Mq of 
hydrogen and type II with more than 0.5 M© of hydrogen. Furthermore we distinguish different cases based on the response of the companion 
star. Case 1: The secondary evolves to an O star at time of explosion of the primary. Case 2: the secondary evolves to a B supergiant. Case 3: the 
secondary evolves to a K supergiant. See also Fig. 4 for an enlargement. 




1000 1200 1400 1600 1800 

Orbital period (days) 

Fig. 4. Zoom-in of Fig. 3 for initial mass ratios close to one. The region (case 3) where the secondary evolves to a K supergiant is subdivided into 
case 3a, where the primary explodes before the secondary, and case 3b, where the secondary explodes before the primary. The location of the four 
models in Table 2) is also indicated. 



of material lost by the primary star as a result of Roche-lobe 
overflow. 

The effect on the primary star is rather modest as its evolu- 
tion has been determined largely by its development before the 
onset of Roche-lobe overflow. In Table 4 we list various prop- 



erties of models computed with accretion efficiencies varying 
between 1 and 0.25. The luminosity and temperature of the pri- 
mary star at the moment of explosion are hardly affected; in all 
cases we find that the primary star will be a cool supergiant at the 
time of explosion. Nevertheless its final mass and therefore the 
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Table 4. Properties of binary stellar models for an initial primary mass of 15 Mq, with an initial orbital period of 1500 days and with variation of 
the initial secondary mass M2i and the accretion efficiency of the companion f}. 





M,,- 
(Mo) 


(Me) 


(mI) 


P 


Mm 
{Me) 




(lO'K) 




(lO^K) 


#1 


14 


5.93 


20.92 


1.00 


0.354 


5.07 


3.48 


4.95 


31.4 


#la 




6.38 


17.23 


0.50 


0.598 


5.07 


3.40 


4.43 


27.1 


#2 


14.35 


5.98 


21.23 


1.00 


0.370 


5.07 


3.47 


5.06 


17.9 


#2a 




6.06 


20.46 


0.90 


0.416 


5.07 


3.45 


5.04 


16.6 


#2b 




6.14 


19.72 


0.80 


0.462 


5.07 


3.43 


5.01 


14.3 


#2c 




6.19 


19.35 


0.75 


0.489 


5.07 


3.43 


4.97 


8.67 


#2d 




6.24 


18.97 


0.70 


0.518 


5.07 


3.42 


4.84 


3.77 


#2e 




6.34 


18.25 


0.60 


0.575 


5.07 


3.40 


4.81 


3.78 


#2f 




6.45 


17.54 


0.50 


0.633 


5.07 


3.39 


4.81 


3.77 


#2g 




6.74 


15.87 


0.25 


0.806 


5.07 


3.37 


4.80 


3.72 


#3 


14.55 


6.01 


21.40 


1.00 


0.390 


5.07 


3.45 


4.85 


3.82 


#3a 




6.49 


17.72 


0.50 


0.660 


5.07 


3.39 


4.85 


3.79 


#A 


8.6 


5.41 


10.44 


0.25 


0.091 


5.07 


3.79 


3.93 


25.4 


#B 


8.8 


5.41 


12.49 


0.50 


0.096 


5.07 


3.80 


4.17 


27.8 


#C 


9.7 


CONTACT 


/ 


0.75 


/ 


/ 


/ 


/ 


/ 


#D 


10 


5.50 


15.49 


0.75 


0.131 


5.07 


3.69 


4.46 


30.5 



Notes. The final mass of the primary and secondary, Mfi and Mf2, and the amount of hydrogen in the envelope Mh are given at the time of explosion, 
together with the luminosity log(L/Lo) and effective temperature Teff of the primary and secondary at the onset of central carbon burning. If the 
evolution of the binary system ended because of formation of a contact binary, the parameters cannot be determined and are indicated by '/'. 
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Fig. 5. Evolution tracks in the Hertzsprung-Russel diagram of the sec- 
ondary of binary systems with initial mass 15 + 14.35 Mg and initial 
orbital period 1500 days. The differences between the evolution tracks 
of the secondary are caused by variation of the accretion efficiency (fi). 
(see models #2 - #2g in Table 4) 

amount of hydrogen in the envelope at the moment of explosion 
Mm increases if we assume a lower accretion efficiency. This 
can be understood as an effect related to response of the orbit 
and thus the size of Roche lobe to mass transfer and mass and 
angular-momentum loss from the system. In our models mass 
loss from the system widens the orbit and therefore the Roche 
lobe. This results in smaller mass-transfer rates and therefore 
larger envelope masses at the time of explosion. This influences 
the parameter space for which we predict type lib supemovae, as 
depicted in Fig. 3. For an accretion efficiency of 50% the range 
of initial orbital periods resulting in SN lib progenitors shifts by 
about 200 days to smaller periods in comparison with the con- 
servative case. 

The evolution and properties of the secondary are more sen- 
sitive to the adopted accretion efficiency. The lower accretion 
rate reduces the expansion of the companion star As a result. 



the formation of contact can be avoided in systems with more ex- 
treme initial mass ratios. This effect widens the parameter space 
for type lib supemovae. In table 4 we list a few test models that 
explore the effect of nonconservative mass transfer for extreme 
mass ratios. For example, we find that the mass ratio leading to 
contact for a system with an initial orbital period of 1500 days 
shifts from Mt/Mi x 0.74 assuming conservative mass transfer 
to about 0.65 if we assume a accretion efficiency of 75% (see 
models #C-#D in Table 4). For lower accretion efficiencies we 
find that formation of contact is no longer the effect that lim- 
its the parameter space for type lib progenitors. In systems with 
such extreme mass ratios, that mass stripping from the primary 
star becomes so efficient that these stars do not have enough hy- 
drogen left at the moment of explosion to become lib. The most 
extreme mass ratio resulting in a lib progenitor assuming the 
same initial orbital period shifts to 0.59 (0.58) assuming an ac- 
cretion efficiency of 50% (25%), see model #B, (#A) in Table 4. 
In section 6 we discuss the influence of y6 on the rate of type lib 
supernova. 

The appearance of the companion star at the moment of 
explosion is also affected. In the first case of accretion during 
the main sequence the secondary rejuvenates, regardless of the 
adopted accretion efficiency. However, the effective temperature 
and luminosity decrease for lower accretion efficiencies as a re- 
sult of the smaller amount of accreted mass (e.g. compare mod- 
els #1 and #la in Table 4). The evolution of the secondary in the 
second case of accretion during the Hertzsprung gap depends 
most strongly on the accretion efficiency. In Fig. 5 we illus- 
trate the evolutionary tracks of the secondary in the Hertspmng- 
Russel diagram. For conservative mass transfer the secondary 
evolves to an overluminous B star The effective temperature de- 
creases with decreasing accretion efficiency. If the accretion ef- 
ficiency is 75%, the secondary evolves to an A supergiant, and 
for accretion efficiencies of 50% or less the secondary evolves to 
a much cooler K supergiant (e.g. compare models #2 and #2g in 
Table 4). In the third case of mass transfer while the secondary 
is already on the giant branch, the secondary still evolves to a K 
supergiant. In this situation its luminosity and effective temper- 
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ature are almost the same as for conservative mass transfer (e.g. 
compare models #3 and #3a in Table 4). 

In summary we conclude that lower accretion efficiencies 
shift and widen the parameter space for type lib supernovae from 
stable mass transfer. However, it proves to be even more difficult 
to produce a companion that is a B supergiant at the moment 
of explosion for non-conservative mass transfer. If the accretion 
efficiency is smaller than about 60% our models predict no com- 
panions to reside in the middle of the Hertzsprung-Russel dia- 
gram. They will either be hot and compact residing in the main 
sequence band or, if the initial masses were very similar, the sec- 
ondary will appear as a K supergiant. 

6. Predicted rates 

The currently most likely observed rate of type lib supemovae 
with respect to all core collapse supernova is about 3 % (see 
Sect. 2 for a discussion). The accuracy of the observed rate is 
expected to go up in the near future thanks to the automated de- 
tection and classification of supernova light curves. Computing 
the rate predicted from our models is not straightforward and re- 
quires adopting further uncertain assumptions about for ex ample 
the binary fraction and the distribution functions of the initial pa- 
rameters. That said, we still consider it worthwhile to estimate 
the rate predicted from our models, which we do below adopt- 
ing commonly made, most reasonable assumptions and assess- 
ing the effect of uncertainties in these assumptions on the rates. 

In the previous sections we discussed progenitor models as- 
suming an initial mass of 15 Mq for the primary star. Based on 
a few test models we conclude that the parameter space does not 
significantly change for systems with primary masses between 
roughly 10 and 20 solar masses. This mass range dominates the 
population of massive binaries in which the primary is massive 
enough to undergo a core collapse supernova. In the following 
derivations we will assume that the characteristics of our conser- 
vative models with a 15 M© primary are representative for the 
population of massive binaries. 

We assume that the initial orbital periods, P, are distributed 
according to Opik's (1924) law (Kouwenhoven et al. 2007), 

/(P)ocp-i for 0.5<P(d)<10^ 

The upper limit is chosen close to the widest orbital period for 
which we expect binaries to interact via Roche-lobe overflow 
during their life time. For the initial mass ratio, which we define 
as the mass of the initially less massive star over the mass of the 
initially most massive star, q s M2/M1, we adopt a power-law 
distribution, 

f{q)ocq' for 0.25 < ^ < 1. 

Binaries with extreme mass ratio are hard to detect and even 
if they have been detected it is still difficult to determine their 
mass ratio accurately. Therefore we limit ourselves to systems 
with M2/M1 > 0.25, following the approach of for example Pols 
et al. (1991). The parameter x describes whether the distribution 
of initial mass ratios is flat {x - 0), skewed towards systems 
equal masses {x > 0) or favors systems with unequal masses 
(x < 0). Although these distribution functions are uncertain, they 
are consistent with observed distributions. Kouwenhoven et al. 
(2007) and references therein find x = -0.4 for the nearby as- 
sociation Scorpius OB2 while x = -1 corresponds to a study of 
spectroscopic binaries by Trimble (1990). A flat or uniform mass 
ratio distribution x - Q has been quoted by Sana et al. (2009) for 
young open clusters, whereas e.g. Pinsonneault & Stanek (2006) 



and Kobulnicky & Fryer (2007) claim that massive binaries like 
to be twins, i.e. x > 0. We will adopt jc = as our standard 
assumption and consider the two extreme cases x = + 1 as well. 

We adopt a binary fraction /bin of 50% as our standard as- 
sumption, i.e. for every single star with mass Mi there is exactly 
one binary system of which the primary mass is Mi and the or- 
bital period and mass ratio are within the ranges specified above. 
For comparison. Mason et al. (1998) derives a binary fraction 
of 59%-75% for O stars in clusters and associations. Garcia 
& Mermilliod (2001) find fractions between 14-80%, whereas 
Kobulnicky & Fryer (2007) infer a fraction higher than 70%. On 
the other hand, Pols et al (1991) adopted a fraction of 27.5% of 
binaries with B-type primaries, P < 10 yr and 0.25 < q' < 1. 
Besides our standard assumption we will consider binary frac- 
tions of /bin = 25% and 100%. We note that exploring the effect 
of different binary fractions partially covers the uncertainties in 
the assumed upper limits for the orbital period and mass ratio 
that we assumed above for the initial parameter distributions. 

The observed rate of type lib supernovae is expressed with 
respect to the rate of all types of core collapse supernova. To 
compute this number we must consider that a significant frac- 
tion of the companion stars end their lives as core collapse su- 
pernovae. The precise fraction will depend on the initial mass 
ratio distribution, the amount of mass lost from the secondary 
and the fraction that is actually accreted by the secondary. For 
clarity, we will just consider the two extreme cases in which all 
or none of the companion stars result in a core collapse super- 
nova. 

6. 1 . The rate of type lib supernovae 

Under the standard assumptions described above, we find that 
the number of type lib supernovae over the number of core col- 
lapse supernovae formed via stable mass transfer as predicted 
from our models is roughly 0.6 %, see Table 5, about a factor 
five lower than the observed rate (Sect. 2). The rate is larger, 
but only by a factor two, than the rate predicted from single 
stars, 0.3%, under the assumption that all stars are single (see 
Section 4). Our predicited rate is compatible with the rate found 
by Podsiadlowski et al. (1992) for the 'stripped' supernovae. 
In contrast to what some previous authors have suggested (e.g. 
Podsiadlowski et al. 1993), we emphasize that not only for the 
single star channel but also for the binary channel one needs to 
fine-tune the initial parameters of the system. For single stars the 
initial mass must be finely-tuned, whereas for binaries the com- 
bination of initial orbital period and inital mass ratio must by be 
carefully arranged. 

The derived rate depends on the adopted initial mass ratio 
distribution, the assumed binary fraction and on whether or not 
the companion star results in a core-collapse supemovae, see 
Table 5. Even when we adopt the assumptions favoring a high 
SN lib rate from binaries (100% binaries, mass ratio distribu- 
tion skewed to equal-mass systems and no core-collapse super- 
novae from the companions) the rate derived from our conserva- 
tive models is still below the observed rate. 

Another uncertainty to consider is the efficiency of mass 
transfer. In our models, lower efficiencies shift and widen the 
parameter space. In Table 5 we show that the number of SNe 
lib over the number of core-collapse supernovae may increase 
to over 1 %, with the warning that for a proper evaluation of the 
parameter space one needs a more extensive model grid. 

A further uncertainty we have to consider is in the range of 
hydrogen-envelope masses in the progenitor that give rise to a 
SN lib. The lower limit of 0. 1 M© is determined by the boundary 
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Table 5. Predicted rates for type lib supemovae from stable mass transfer in binaries and its dependence on different input assumptions. 



Mass ratio distr. fav. extreme ratios 


flat 


fav. equal 


masses 


(x = -l) 


(x = 0) 


(x = 


1) 



Fraction of binaries producing type lib SNe via stable mass transfer 

0.83% 1.33% 1.85% 

Number of type lib relative to the number of all core collapse SNe 
(assuming that the companion does - does not produce a core collapse SNe) 

/b,„ = 25% 0.17 - 0.21% 0.27 - 0.33% 0.37 - 0.46% 

fb,„ = 50% 0.28-0.42% 0.44-0.67% 0.62-0.93% 
/bin = 100% 0.42 - 0.83% 0.67 - 1.33% 0.93 - 1.85% 

Estimates of the rates assuming lower accretion efficiencies /?' 

/? = 1.0 0.44-0.67% 
P = 0.75 0.60 - 0.92 Vc^ 

13 = 0.5 0.69- 1.06%' 

/3 = 0.25 0.71 - 1.09%^ 

Assuming that all systems with orbital periods between 1000 and 2000 days produce lib's'-' 

2.3-3.5% 



'Assuming a binary fraction of /bi„ = 50%. 

^The estimates for lower accretion efficiencies are based on just a few test models and should only be taken as indicative. A more extended grid of 

non-conservative models is required for a proper evaluation of these rates. 

'This estimate does no longer depend on the assumed distribution of initial mass ratios. 

Notes. See Section 6 for details. The rate predicted based on our standard assumptions is marked in boldface. 

Table 6. Relative rates for the properties of the companion star at the moment of explosion for different initial mass ratio distributions based on 
conservative models. 



Mass ratio distr. 


fav. extreme ratios flat 


fav. equal masses 




(x = -l) (x = 0) 


(x= 1) 


Case 1 (O star) 


91.5% 90.3% 


89.1% 


Case 2 (B supergiant) 


2.5% 2.8% 


3.2% 


Case 3 (K supergiant) 


6.0% 6.9% 


7.7% 



between compact and extended type lib. The upper limit of about 
0.5 M0 is determined by the lack of a plateau in the light curve, 
which depends on the chemical structure of the stellar envelope 
and therefore is less certain. If we change the boundary from 0.5 
to 0.6 M0, test models show that the initial orbital period range 
will only widen by about 80 days. This means an increase of 
the overal rate of type lib SNe of about 10%. The uncertainties 
in these boundaries will therefore not have a large effect on the 
rate. 

As an extreme assumption, we consider that all systems with 
orbital periods between 1000 and 2000 days produce lib pro- 
genitors. This corresponds to 7% of all binary systems with ini- 
tial mass ratios and orbital periods in the ranges specified ear- 
lier. This would increase the relative rate compared to the rate of 
core collapse supernova to 2.3-3.5%, depending on whether the 
companion star produces a core-collapse supernova or not and 
assuming a binary fraction of 50%. This number is consistent 
with the observed rate. 



6.2. The relative rates of different cases for tine companion 

Table 6 indicates the relative rates for the different characteris- 
tics of the companion at the moment of explosion, as described 
in section 5.1. The large majority, about 90%, of type lib su- 
pernova resulting from stable mass transfer is expected to have 
an O-star companion at the moment of explosion. The rate of 



systems with a blue supergiant companion, which applies to SN 
1993 J and possibly 200 lig, is predicted to be quite low accord- 
ing to ourmodels, about 3% of the type lib supemovarate. These 
relative rates reflect the large and very small regions in parameter 
space for O-star and B-supergiant companions, respectively, as 
depicted in Fig. 3. However, these percentages do not necessar- 
ily reflect the probability of detecting a companion of a certain 
type. Detecting an O-star companion in post-explosion images 
of a supernova is more difficult than detecting a B or K super- 
giant, as we demonstrated in Section 5.1. 

The probability for the presence of blue companions to 
type lib SNe we derive is strikingly lower than predicted by 
Podsiadlowski et al. (1992). They find that, if the companion 
starts to accrete after becoming a giant star, it will contract in 
a similar way to a feature known as "blue loops" seen in some 
evolutionary tracks of single stars. We do not find this behavior 
with our stellar evolution code. Whether or not stellar models 
perform blue loops is very sensitive to the ratio of core mass to 
total stellar mass and to details in the chemical profile outside the 
stellar core (Kippenhahn & Weigert 1994). If a significant frac- 
tion of the case 3 models would result in blue companions, the 
relative rate for the presence of blue companions may increase 
by a factor of two or three (see Table 6). 

In the most common case, accretion during the main se- 
quence, the star rejuvenates, i.e. the size of its convective core 
increases to adapt to the new stellar mass, see section 5.1. In our 



11 



J.S.W. Claeys et al.: Binary progenitor models of type lib supemovae 



code, all main sequence stars that accrete rejuvenate. However, 
whether stars rejuvenate or not depends on the assumed effi- 
ciency of semi-convection. Stars that have evolved towards the 
end of the main sequence build up a chemical gradient between 
the helium rich core and the hydrogen rich envelope. Braun & 
Langer (1995) show that, assuming a low efficiency of semi- 
convective mixing, the growth of the convective core in an ac- 
creting main-sequence star is prevented. These stars may also 
appear as blue supergiants at the moment of explosion. 

Based on their models, we estimate that in our case of con- 
servative mass transfer the companion must have reached a cen- 
tral helium mass fraction of about 0.75 or more to prevent reju- 
venation. In our models this occurs for initial mass ratios larger 
than about 0.87, and we estimate that this may increase the rel- 
ative rate of the presence of blue companions up to about 50%. 
For lower accretion efficiencies it becomes even easier to prevent 
rejuvenation: a central helium mass fraction for the accreting star 
of just 0.6 and higher are required if ;6 = 0. However, the com- 
panion must accrete enough to finish its main sequence evolu- 
tion before the explosion of the primary star, to be observed as a 
blue supergiant. More detailed estimates are beyond the scope of 
this study, but we do emphasize this as an interesting possibility 
to use the companions of supernovae to gain insight in internal 
mixing processes. 

7. Discussion & Conclusion 

We identified binary progenitor models for extended type lib su- 
pernovae. In these models the most massive star is stripped from 
its envelope by interaction with its companion, such that only 
several tenths of solar masses of hydrogen remain at the mo- 
ment of explosion. We find that the most massive star must fill 
its Roche lobe during central helium burning in order to achieve 
this. We derive for which range of initial orbital periods and ini- 
tial mass ratios binary systems are expected to produce type lib 
supernova progenitors. 

We have discussed in detail the properties of the companion 
star at the moment of explosion, motivated by the detection of a 
companion for SN 1993 J, a possible companion for 200 lig and 
the non-detection of a companion in the supernova remnant of 
Cas A. We distinguish three cases: (1) the companion appears 
as a hot O star, when it accreted while still being on the main 
sequence, (2) the companion becomes an over-luminous B-star 
if mass transfer started while the companion was crossing the 
Hertzsprung gap, a scenario which may apply to 1993 J and pos- 
sibly 2001ig and (3) the companion will be a K-supergiant when 
it accreted after ascending the giant branch. The third case ap- 
plies to systems with very similar initial masses and can be sub- 
divided in (3a) systems in which the initially most massive star 
explodes first and (3b) systems in which the companion star ex- 
plodes first, up to 10^ years before the primary. These models 
may apply to Cassiopeia A -if it was an extended-lib SN- ex- 
plaining the lack of evidence for a companion^. If the progenitor 
of Cas A was compact, not a red supergiant, at time of explosion 
another evolutionary scenario then discussed here could explain 
this specific supernova remnant, e.g. a common envelope sce- 
nario. 

However, our models predict that the scenarios we propose 
for 1993 J, 2000ig (case 2) and Cas A (case 3b) are very rare and 

- Signatures from the first supernova, for example in radio emission 
or from interaction with the interstellar medium will be hard to detect as 
Cas A is bright radio source itself and the remnant of the first supernova 
will be roughly 10 to 100 times larger than the remnant of Cas A. 



require that the companion efficiently accretes the transferred 
material. The accretion efficiency, i.e. the fraction of transferred 
material that is accreted by the companion, is a major uncer- 
tainty for binary evolutionary models. A lower accretion effi- 
ciency shifts the parameter space to smaller orbital periods (by 
about 200 days for an efficiency of 50% compared to conser- 
vative mass transfer) and widens the parameter space towards 
more extreme initial mass ratios (assuming an efficiency of 50% 
we estimate that the lib rate predicted by our models increases 
by roughly a factor 1.6 compared to conservative mass transfer). 
However a B supergiant is identified more easily in compari- 
son with the other possible evolutionary paths of the companion 
star, as pointed out in section 5.1 and Table 3, which gives a bias 
towards a higher fraction of observed B-supergiants as a com- 
panion star 

Our conservative models predict ~ 0.6% of all core collapse 
supernova to be of type lib (using our standard assumptions, 
see Section 6), about a factor of five lower than the currently 
most likely observed rate. This rate is larger, but only by a factor 
of two than the rate predicted for a pure single star population, 
~ 0.3%, in which high mass stars lose their envelope due to stel- 
lar winds. We emphasize that both the binary and single stars 
scenario only produce type lib in a very limited region of the 
initial parameter space. 

We warn the reader that the rates quoted above are uncertain. 
The rate from single stars is affected by our limited understand- 
ing of mass loss rates from massive red giants. The rate from 
binary stars is affected by uncertainties in the physics of binary 
interaction as well as the adopted distributions of initial binary 
parameters. If we make several assumptions that favor a high 
rate of lib supernovae (a high binary fraction, a distribution of 
initial mass ratios which is skewed towards systems with equal 
masses, a low mass-transfer efficiency and only moderate angu- 
lar momentum loss from the system), we derive a rate that is 
consistent with the lower limit of the currently most reliable ob- 
served rate of Smartt et al. (2009). 

In addition, the observed rate is based on only a few SNe and 
it is difficult to make a distinction between an extended type lib 
and a compact, the latter of which is more similar to a type lb 
SN and is formed by another evolutionary scenario (e.g. Yoon 
et al. 2010) . Therefore, the observed rate is also one of the un- 
certainties. Nevertheless, the rates predicted from our standard 
models seem to indicate that there is definitely room for a single 
star scenario and a common envelope scenario leading to a type 
lib. 

A brief comparison witti previous work 

Our findings are consistent with earlier studies by Podsiadlowski 
et al. (1992, 1993), Woosley et al. (1994), Maund et al. (2004) 
and Stancliffe & Eldridge (2009). Nevertheless we can note two 
major differences. Contrary to suggestions by Podsiadlowski 
et al. (1993) we find that it proves to be hard to explain the 
presence of a blue companion star at the time of explosion, con- 
firming the thought raised by Stancliffe & Eldridge (2009) with 
our more extended model grid. In addition, we find that the pri- 
mary star needs to fill its Roche lobe during helium burning (late 
Case B mass transfer), whereas Podsiadlowski et al. (1993) re- 
quire mass transfer to start after helium burning (Case C mass 
transfer). This can be attributed to differences in the adopted 
physics: the adopted initial composition, the opacity tables and 
the mixing length parameter which have a not negligible effect 
on the radii of giants and therefore on the onset of and response 
to Roche-lobe overflow. Furthermore there are small differences 
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in prescription for the mass-transfer rate, the treatment of non- 
conservative mass transfer and the mass loss occurring between 
the last computed model and the explosion. As a result we find 
that slightly different initial orbital periods are required to pro- 
duce a lib progenitor However, if we allow for uncertainties, our 
estimate of the lib rate is consistent with previous studies. 

What type lib SNe teach us about stellar and binary physics 

The rate of type lib supemovae and the relative occurrence of 
different companions have the potential to constrain both stel- 
lar physics, such as internal mixing processes and stellar wind 
mass loss, but also the physics of interacting binaries, when the 
accuracy of the determined rates increases in the near future. 

The discrepancy between the currently most likely observed 
rate and the rate predicted by our models seems to point towards 
lower accretion efficiencies. At lower accretion efficiencies, sys- 
tems with more extreme mass ratios can avoid contact and re- 
tain enough hydrogen in their envelope to become progenitors 
of lib supemovae. On the other hand, the presence of a blue 
companion for SN 1993 J and possibly SN 2000ig indicates that 
at least a substantial fraction must be accreted. Therefore, im- 
proving statistics about lib supemovae may in the future help to 
constrain the efficiency of mass transfer, one of the major uncer- 
tainties in binary evolutionary models. 

The need to avoid contact and a subsequent spiral in for these 
binaries with an evolved giant companion also gives information 
about two other uncertain input assumption related to the physics 
of mass transfer It may indicate that the entropy of accreted ma- 
terial must be low, to prevent swelling of the companion, and 
it may indicate that the mass-transfer rate in such systems re- 
mains lower than the self-regulated thermal-timescale rate, see 
section 3. 

Finally the detection of one and possibly two blue compan- 
ions while our standard models predict that this should be a rare 
event is puzzling. As discussed in section 6.2, this may be ex- 
plained with a low efficiency of semi-convection. 

Outlook 

We illustrated the potential of type lib supemovae, being on the 
border between the hydrogen rich type 11 and the hydrogen poor 
type Ib/c supernova, to give insight in stellar physics such as 
internal mixing processes and stellar winds as well as physics 
of interacting binaries. Even though the predictions from stellar 
evolutionary models are still challenged by various uncertainties, 
it has become feasible to compute large grids of binary evolu- 
tionary models and assess the impact of the various assumptions. 
Especially from the observational side we expect large progress 
in the near future. Statistics are improving with the large sam- 
ples resulting from automated surveys. These will allow us to 
even asses the rate of lib supemovae in different host galaxies 
probing different environments and even metallicity regimes. 
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Appendix A: Tables 



Table A.l. Properties of binary stellar models for an initial primary mass of 15 Mq, conservative mass transfer and with variation of the initial 
secondary mass M2i and initial orbital period Po,b. 
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1 1 nn 
1 lUU 


O.izo/ 


C AC 


Q AA 

3.69 


C 1 c 

J. 15 


13.3 


id 


1 A on 
14.6 / 


5.31 


ZZ. ol 


1 1 nn 
1 lUU 


O.izoz 


C AC 


Q AA 

3.69 


C 1 c 

5.15 


1 A A 

14.4 


le 


1 A 

14.66 


5.31 


ZZ. oO 


1 1 nn 
1 lUU 


0.1 zoO 


C AC 


Q AA 

3.09 


C 1 c 

5.15 


1 C A 

15.9 


za 


1 A m 
14.93 


5.4o 


ZZ.4Z 


1 onn 
IZUU 


0.15/0 


C AA 

5. Do 


Q AA 
3. DO 


C 1 c 

5.15 


1 Q y1 
13.4 


ZD 


1 A m 
14.92 


5.4d 


ZZ. 41 


1 onn 
IzUU 


0.14/0 


C AA 

5.U6 


Q AA 

3.66 


C 1 c 

5.15 


1 Q A 

13.6 


zc 


1 /I o 

14. y 


J. 40 


TO /in 
ZZ.4U 


1 onn 
IZUU 


n 1 CCA 
U.l JJO 


C AA 
J.UO 


1 AA 
3. DO 


C 1 c 

5.15 


13.3 


zd 


1 A oc 

14.65 


Z AZ 

5.45 


ZZ.35 


1 onn 
IZUU 


n 1 cc c 
0.15jj 


C AA 

5. Do 


Q AT 

3.0 / 


C 1 /I 

5. 14 


13. / 


ze 


1 A OA 

14.64 


Z AZ 

5.45 


ZZ.34 


1 onn 
IZUU 


n 1 cc T 
0.1553 


C AA 

5. Do 


Q AT 

3.0 / 


C 1 /I 

5. 14 


1 C A 

15.4 


3a 


14.73 


5.62 


22.00 


1300 


0.1949 


5.07 


3.60 


5.12 


16.3 


3d 


1 A '70 

14. /z 


5.0Z 


ZZ.OO 


1 inn 
13UU 


n 1 nco 
0.1956 


C A"7 

5.0/ 


Q AA 

3.DU 


C 1 1 

5.12 


I A A 

16.0 


3c 


14.71 


5.61 


21.99 


1300 


0.191 1 


5.07 


3.60 


5.12 


14.9 


3d 


14.7 


5.61 


21.98 


1300 


0.191 1 


5.07 


3.60 


5.12 


14.9 


3e 


1 A /^O 

14.00 


5.61 


Zl.96 


1 inn 
13UU 


0.1912 


C A"? 

5.0/ 


Q AA 
3.DU 


C 1 1 

5.11 


ICO 

15.2 




14 (S6 


J.U 1 


91 Q4 
z 1 .y^ 




U.^O.iO 


■S 07 


J.UO 


5.10 


1 8 4 


4a 


14.54 


5.79 


21.63 


1400 


0.2774 


5.07 


3.53 


5.09 


17.6 


4b 


14.53 


5.79 


21.62 


1400 


0.2787 


5.07 


3.53 


5.09 


16.6 


4c 


14.5 


5.79 


21.59 


1400 


0.2711 


5.07 


3.53 


5.08 


17.1 


5a 


14.36 


5.98 


21.23 


1500 


0.3677 


5.07 


3.47 


5.06 


18.1 


5b 


14.35 


5.98 


21.23 


1500 


0.3699 


5.07 


3.47 


5.06 


17.9 


5c 


14.33 


5.98 


21.21 


1500 


0.3625 


5.07 


3.47 


5.06 


18.2 


5d 


14.32 


5.97 


21.20 


1500 


0.3765 


5.07 


3.47 


5.04 


19.7 


6a 


14.17 


6.20 


20.80 


1600 


0.4772 


5.07 


3.42 


5.03 


18.8 


6b 


14.15 


6.19 


20.79 


1600 


0.4749 


5.07 


3.43 


5.02 


19.1 


7a 


14.01 


6.44 


20.37 


1700 


0.6257 


5.07 


3.39 


5.00 


19.1 


7b 


14 


6.44 


20.37 


1700 


0.6458 


5.07 


3.39 


5.00 


19.0 


7c 


13.97 


6.43 


20.34 


1700 


0.6399 


5.07 


3.39 


4.99 


20.0 



Notes. The secondary is, for every model, at the start of mass transfer in its Hertszprung gap. The final mass of the primary and secondary, Mfi 
and Mf2, and the amount of hydrogen in the envelope Mh are given at the time of explosion, together with the luminosity log(L/LQ) and effective 
temperature Teff of the primary and secondary at the onset of central carbon burning. 
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Table A.2. Properties of binary stellar models for an initial primary mass of 15 Mq, conservative mass transfer and with variation of the initial 
secondary mass Mai and initial orbital period Pet. 





(Mo) 


(Mo) 


My 

(Mo) 


Port 

(Days) 


Mm 
(Mo) 




(lO'K) 




7'eff,2 

(lO^K) 


hOa 


13.5 


5.10 


21.52 


1000 


0.0835 


5.03 


3.86 


4.96 


32.2 


hOb 


14 


5.13 


21.98 


1000 


0.0886 


5.03 


3.82 


5.01 


31.7 


hOc 


14.85 


5.19 


22.76 


1000 


0.0980 


5.04 


3.76 


5.09 


30.3 


hla 


13.3 


5.16 


21.25 


1100 


0.0997 


5.03 


3.76 


4.94 


32.3 


hlb 


13.4 


5.17 


21.34 


1100 


0.1014 


5.04 


3.76 


4.95 


32.2 


hlc 


14 


5.23 


21.85 


1100 


0.1111 


5.04 


3.73 


5.00 


31.7 


hid 


14.85 


5.31 


22.59 


1100 


0.1258 


5.05 


3.69 


5.08 


30.2 


h2a 


13.2 


5.31 


20.86 


1200 


0.0983 


5.05 


3.78 


4.92 


32.2 


h2b 


13.3 


5.32 


20.95 


1200 


0.1012 


5.04 


3.77 


4.93 


32.1 


h2c 


14 


5.37 


21.60 


1200 


0.1195 


5.06 


3.72 


4.99 


31.5 


h2d 


14.8 


5.45 


22.31 


1200 


0.1450 


5.06 


3.67 


5.07 


30.1 


h2e 


14.83 


5.45 


22.33 


1200 


0.1450 


5.06 


3.67 


5.07 


30.1 


h3a 


12.1 


5.35 


19.65 


1300 


0.0983 


5.06 


3.78 


4.80 


32.6 


h3b 


12.2 


5.36 


19.75 


1300 


0.1079 


5.06 


3.77 


4.81 


32.5 


h3c 


14.64 


5.61 


21.93 


1300 


0.2009 


5.07 


3.60 


5.04 


30.5 


h4a 


11.4 


CONTACT 


/ 


1400 


/ 


/ 


/ 


/ 


/ 


h4b 


11.5 


5.43 


18.85 


1400 


0.1223 


5.06 


3.72 


4.74 


32.5 


h4c 


11.9 


5.47 


19.29 


1400 


0.1329 


5.06 


3.69 


4.77 


32.6 


h4d 


12 


5.48 


19.38 


1400 


0.1430 


5.07 


3.68 


4.78 


32.6 


h4e 


14.49 


5.79 


21.58 


1400 


0.2785 


5.07 


3.53 


5.01 


30.8 


h5a 


11.1 


CONTACT 


/ 


1500 


/ 


/ 


/ 


/ 


/ 


h5b 


11.2 


5.56 


18.37 


1500 


0.1615 


5.07 


3.65 


4.70 


32.4 


h5c 


11.6 


5.60 


18.84 


1500 


0.1843 


5.07 


3.62 


4.73 


32.5 


h5d 


13 


5.79 


20.06 


1500 


0.2711 


5.07 


3.53 


4.86 


32.3 


h5e 


14 


5.93 


20.92 


1500 


0.3538 


5.07 


3.48 


4.95 


31.4 


h5f 


14.31 


5.97 


21.19 


1500 


0.3604 


5.07 


3.47 


4.98 


31.0 


h6a 


10.7 


CONTACT 


/ 


1600 


/ 


/ 


/ 


/ 


/ 


h6b 


10.8 


5.68 


17.81 


1600 


0.2183 


5.07 


3.58 


4.64 


32.4 


h6c 


11.3 


5.76 


18.36 


1600 


0.2442 


5.07 


3.55 


4.69 


32.5 


h6d 


13 


6.01 


19.81 


1600 


0.3916 


5.07 


3.46 


4.84 


32.2 


h6e 


14.1 


6.19 


20.73 


1600 


0.4740 


5.07 


3.43 


4.94 


31.2 


h6f 


14.14 


6.19 


20.77 


1600 


0.4862 


5.07 


3.43 


4.95 


31.1 


h7a 


10.4 


CONTACT 


/ 


1700 


/ 


/ 


/ 


/ 


/ 


h7b 


10.5 


5.81 


17.34 


1700 


0.2716 


5.08 


3.53 


4.61 


32.1 


h7c 


11.1 


5.91 


17.98 


1700 


0.3318 


5.08 


3.50 


4.66 


32.3 


h7d 


11.4 


5.96 


18.23 


1700 


0.3468 


5.08 


3.48 


4.68 


32.4 


h7e 


12.7 


6.20 


19.31 


1700 


0.4784 


5.07 


3.42 


4.80 


32.2 


h7f 


12.8 


6.21 


19.38 


1700 


0.5041 


5.07 


3.42 


4.81 


32.2 


h8a 


10.3 


CONTACT 


/ 


1800 


/ 


/ 


/ 


/ 


/ 


h8b 


10.4 


5.95 


17.08 


1800 


0.3467 


5.08 


3.49 


4.58 


32.0 


h8c 


11 


6.05 


17.70 


1800 


0.4006 


5.08 


3.45 


4.63 


32.4 


hSd 


11.5 


6.16 


18.11 


1800 


0.4668 


5.08 


3.43 


4.68 


32.2 


h8e 


11.8 


6.22 


18.34 


1800 


0.4860 


5.08 


3.42 


4.70 


32.2 


h8f 


11.9 


6.24 


18.41 


1800 


0.5104 


5.08 


3.42 


4.71 


32.2 


h9a 


10.3 


CONTACT 


/ 


1900 


/ 


/ 


/ 


/ 


/ 


h9b 


10.4 


6.12 


16.95 


1900 


0.4531 


5.08 


3.44 


4.57 


31.8 


h9c 


10.8 


6.19 


17.24 


1900 


0.4873 


5.08 


3.42 


4.60 


28.8 


h9d 


10.9 


6.22 


17.31 


1900 


0.5014 


5.08 


3.42 


4.61 


32.0 



Notes. The secondary is, for every model, at the start of mass transfer a main sequence star. The other symbols have a similar meaning as in Table 
A.l. If the evolution of the binary system ended because of formation of a contact binary, the parameters cannot be determined and are indicated 
by 7. 
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Table A.3. Properties of binary stellar models for an initial primary mass of 15 Mq, conservative mass transfer and with variation of the initial 
secondary mass Mai and initial orbital period Pcb. 





M2, 
(Mo) 


Ml,- 
(Mo) 


M2f 

(Mo) 


Port 

(Days) 


M„i 
(Mo) 




(lO'K) 




T<:S.l 

(lO'K) 


*1/ * 2 


time 

(yr) 


g2a 


14.94 


5.60 


22.26 


1200 


0.2178 


5.06 


3.66 


5.08 


3.63 


*2 


14200 


g2b 


14.99 


5.61 


21.75 


1200 


0.2211 


5.06 


3.66 


5.08 


3.63 


*2 


70400 


g3a 


14.74 


5.62 


22.01 


1300 


0.2061 


5.07 


3.60 


4.90 


3.79 


*1 


/ 


g3b 


14.75 


5.62 


22.02 


1300 


0.2058 


5.07 


3.60 


4.90 


3.79 


*1 


/ 


g3c 


14.93 


5.64 


22.18 


1300 


0.2036 


5.07 


3.59 


5.02 


3.67 


*1 


/ 


g3d 


14.94 


5.76 


22.18 


1300 


0.2638 


5.07 


3.59 


5.07 


3.64 


*2 


300 


g3e 


14.99 


5.80 


21.37 


1300 


0.2894 


5.07 


3.58 


5.08 


3.62 


*2 


63100 


g4a 


14.55 


5.79 


21.64 


1400 


0.2761 


5.07 


3.53 


4.86 


3.83 


*1 


/ 


g4b 


14.9 


5.84 


21.94 


1400 


0.3081 


5.07 


3.51 


4.95 


3.74 


*1 


/ 


g4c 


14.94 


5.84 


21.98 


1400 


0.3087 


5.07 


3.51 


5.04 


3.66 


*1 


/ 


g4d 


14.95 


6.02 


21.79 


1400 


0.4012 


5.07 


3.51 


5.07 


3.63 


*2 


7900 


g4e 


14.99 


6.03 


21.01 


1400 


0.4012 


5.07 


3.51 


5.08 


3.61 


*2 


58300 


g5a 


14.37 


5.98 


21.24 


1500 


0.3656 


5.07 


3.47 


4.81 


3.85 


*1 


/ 


g5b 


14.4 


5.99 


21.27 


1500 


0.3728 


5.07 


3.47 


4.82 


3.85 


*1 


/ 


g5c 


14.55 


6.01 


21.40 


1500 


0.3904 


5.07 


3.45 


4.85 


3.82 


*1 


/ 


g5d 


14.9 


6.06 


21.69 


1500 


0.4240 


5.07 


3.45 


4.95 


3.74 


*1 


/ 


g5e 


14.94 


6.06 


21.73 


1500 


0.4244 


5.07 


3.45 


5.04 


3.66 


*1 


/ 


g5f 


14.95 


6.27 


21.56 


1500 


0.5305 


5.07 


3.45 


5.07 


3.63 


*2 


9100 


g5g 


14.99 


6.28 


20.56 


1500 


0.5328 


5.07 


3.45 


5.07 


3.60 


*2 


59500 


g6a 


14.18 


6.20 


20.81 


1600 


0.4796 


5.07 


3.42 


4.78 


3.88 


*1 


/ 


g6b 


14.2 


6.20 


20.83 


1600 


0.4980 


5.07 


3.42 


4.77 


3.88 


*1 


/ 


g6c 


14.3 


6.22 


20.91 


1600 


0.4997 


5.07 


3.42 


4.80 


3.86 


*1 


/ 


g6d 


14.4 


6.23 


20.99 


1600 


0.5236 


5.07 


3.42 


4.82 


3.85 


*1 


/ 



Notes. The secondary is, for every model, at the start of mass transfer on the giant branch. The column with *l/*2 indicates which star in the 
binary system explodes first. The time difference between the two stars reaching carbon burning is indicated if star 2 explodes before star 1. The 
other symbols have a similar meaning as in Table A. 1. 
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Table A.4. Properties of binary stellar models for an initial primary mass of 15 Mq, with variation of the initial secondary mass M2i, variation of 
the accretion efficiency [j and initial orbital period Po,b- 





(Mo) 


Mif 
(Me) 


M2f 

(Ma) 


Porb 

(Days) 


P 


Mm 
(Me) 




(lO^K) 




(lO'K) 


*1/ *2 


t time 

(yr) 


8.5 


14.9 


5.15 


18.86 


800 


0.50 


0.0897 


5.04 


3.84 


5.08 


3.56 


*1 


/ 


9.9 


14.9 


5.13 


22.07 


900 


0.90 


0.0814 


5.03 


3.84 


5.14 


12.6 


*1 


/ 


9.5 


14.9 


5.38 


18.75 


900 


0.50 


0.1671 


5.04 


3.71 


5.10 


3.55 


*2 


11900 


0c9 


14.9 


5.23 


21.97 


1000 


0.90 


0.1050 


5.04 


3.73 


5.13 


11.8 


*1 


/ 


lb9 


14.95 


5.50 


21.41 


1100 


0.90 


0.2107 


5.05 


3.66 


5.08 


3.62 


*2 


20200 


lb7 


14.95 


5.60 


19.93 


1100 


0.70 


0.2469 


5.05 


3.62 


5.07 


3.59 


*2 


16700 


lc9 


14.9 


5.36 


21.81 


1100 


0.90 


0.1512 


5.05 


3.66 


5.12 


10.7 


*1 


/ 


lc5 


14.9 


5.73 


18.54 


1100 


0.50 


0.2878 


5.06 


3.58 


5.09 


3.55 


*2 


12300 


ld75 


14.87 


5.43 


20.56 


1100 


0.75 


0.1712 


5.05 


3.63 


5.08 


8.81 


*1 


/ 


ld5 


14.87 


5.57 


18.57 


1100 


0.50 


0.2680 


5.06 


3.58 


4.94 


3.69 


*1 


/ 


ld25 


14.87 


5.74 


16.67 


1100 


0.25 


0.2857 


5.06 


3.53 


4.93 


3.64 


*1 


/ 


2c5 


14.9 


5.76 


18.47 


1200 


0.50 


0.2661 


5.07 


3.54 


4.95 


3.66 


*1 


/ 


3d5 


14.7 


5.96 


18.16 


1300 


0.50 


0.3658 


5.07 


3.48 


4.89 


3.71 


*1 


/ 


4a7 


14.54 


6.01 


19.34 


1400 


0.70 


0.3906 


5.07 


3.46 


4.85 


3.77 


*1 


/ 


4c9 


14.5 


5.96 


19.68 


1400 


0.90 


0.3671 


5.07 


3.47 


5.02 


12.7 


*1 


/ 


4c7 


14.5 


6.00 


19.31 


1400 


0.70 


0.3873 


5.07 


3.46 


5.00 


9.64 


*1 


/ 


4c5 


14.5 


6.18 


17.84 


1400 


0.50 


0.4844 


5.07 


3.43 


4.86 


3.74 


*1 


/ 


g4b7 


14.9 


6.06 


19.66 


1400 


0.70 


0.4204 


5.07 


3.45 


4.95 


3.69 


*1 


/ 


g4c7 


14.94 


6.07 


19.70 


1400 


0.70 


0.4191 


5.07 


3.45 


5.04 


3.61 


*1 


/ 


g4d7 


14.95 


6.28 


19.60 


1400 


0.70 


0.5387 


5.07 


3.45 


5.07 


3.58 


*2 


9700 


h5.0 


8.3 


5.37 


10.15 


1500 


0.25 


0.0795 


5.07 


3.87 


3.88 


25.1 


*1 


/ 


h5.1 


8.6 


5.41 


10.44 


1500 


0.25 


0.0912 


5.07 


3.79 


3.93 


25.4 


*1 


/ 


h5.2 


8.8 


5.44 


10.64 


1500 


0.25 


0.1074 


5.07 


3.75 


3.96 


25.5 


*1 


/ 


h5.3 


8.8 


5.41 


12.49 


1500 


0.50 


0.0961 


5.07 


3.80 


4.17 


27.8 


*1 


/ 


h5.4 


9.1 


5.44 


12.78 


1500 


0.50 


0.1164 


5.07 


3.75 


4.20 


28.0 


*1 


/ 


h5.5 


9.4 


5.55 


11.22 


1500 


0.25 


0.1510 


5.07 


3.65 


4.06 


26.0 


*1 


/ 


h5.6 


9.4 


5.49 


13.07 


1500 


0.50 


0.1239 


5.07 


3.70 


4.24 


28.2 


*1 


/ 


h5.7 


9.4 


CONTACT 


/ 


1500 


0.75 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


h5.8 


9.7 


CONTACT 


/ 


1500 


0.75 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


h5.9 


10 


5.50 


15.49 


1500 


0.75 


0.1308 


5.07 


3.69 


4.46 


30.5 


*1 


/ 




1 A 

14 


6.36 


1 i.Ii 


ijUO 


U.DO 


0.D9/D 


J.U/ 


Q /in 
3.4U 


4.4o 


I /.I 




/ 


5b9 


14.35 


6.06 


20.46 


1500 


0.90 


0.4163 


5.07 


3.45 


5.04 


16.6 


*1 


/ 


5b8 


14.35 


6.14 


19.72 


1500 


0.80 


0.4624 


5.07 


3.43 


5.01 


14.3 


*1 


/ 


5b75 


14.35 


6.19 


19.35 


1500 


0.75 


0.4886 


5.07 


3.43 


4.97 


8.67 


*1 


/ 


5b7 


14.35 


6.24 


18.97 


1500 


0.70 


0.5176 


5.07 


3.42 


4.84 


3.77 


*1 


/ 


5b6 


14.35 


6.34 


18.25 


1500 


0.60 


0.5750 


5.07 


3.40 


4.81 


3.78 


*1 


/ 


5b5 


14.35 


6.45 


17.54 


1500 


0.50 


0.6325 


5.07 


3.39 


4.81 


3.77 


*1 


/ 


5b25 


14.35 


6.74 


15.87 


1500 


0.25 


0.8062 


5.07 


3.37 


4.80 


3.72 


*1 


/ 


g5c5 


14.55 


6.49 


17.72 


1500 


0.50 


0.6595 


5.07 


3.39 


4.85 


3.79 


*1 


/ 


6b9 


14.15 


6.29 


20.04 


1600 


0.90 


0.5445 


5.07 


3.41 


5.00 


17.5 


*1 


/ 



Notes. The column with *l/*2 indicates which star in the binary system explodes first. The time difference between the two stars reaching carbon 
burning is indicated if star 2 explodes before star 1 . If the evolution of the binary system ended because of formation of a contact binary, the 
parameters cannot be determined and are indicated by '/'• The other symbols have a similar meaning as in Table A.l. 



